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ABSTRACT 

Understanding how disks dissipate is essential to studies of planet formation. 
However, identifying exactly how dust and gas dissipates is complicated due to 
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difficulty in finding objects clearly in the transition of losing their surrounding 
material. We use Spitzer IRS spectra to examine 35 photometrically-selected 
candidate cold disks (disks with large inner dust holes). The infrared spectra 
are supplemented with optical spectra to determine stellar and accretion prop- 
erties and 1.3mm photometry to measure disk masses. Based on detailed SED 
modeling, we identify 15 new cold disks. The remaining 20 objects have IRS 
spectra that are consistent with disks without holes, disks that are observed 
close to edge-on, or stars with background emission. Based on these results, we 
determine reliable criteria for identifying disks with inner holes from Spitzer pho- 
tometry and examine criteria already in the literature. Applying these criteria to 
the c2d surveyed star-forming regions gives a frequency of such objects of at least 
4% and most likely of order 12% of the YSO population identified by Spitzer. 

We also examine the properties of these new cold disks in combination with 
cold disks from the literature. Hole sizes in this sample are generally smaller 
than for previously discovered disks and reflect a distribution in better agreement 
with exoplanet radii. We find correlations between hole size and both disk and 
stellar masses. Silicate features, including crystalline features, are present in the 
overwhelming majority of the sample although 10 /zm feature strength above the 
continuum declines for holes with radii larger than ~7 AU. In contrast, PAHs are 
only detected in 2 out of 15 sources. Only a quarter of the cold disk sample shows 
no signs of accretion, making it unlikely that photoevaporation is the dominant 
hole forming process in most cases. 

Subject headings: stars: planetary systems: protoplanetary disks - stars: pre- 
main sequence 



1. Introduction 

Near- and mid-IR observations of young stars demonstrate that optically-thick circum- 



are nearly entirely absent around members of 10 Myr old associations (e.g. . 
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i et al. 


2001: Gutermuth et al. 2004: Sicilia-Aeuilar et al. 2005. Low et al. 


2005; 
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rrie et al. 
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Accretion ceases on approximately the same timescale (e.g. 


Calvet et al. 


200, c 


i). The dis- 



appearance of gas and dust - planet ary building materi a l - places stringent limits on the 



timescales of giant planet formation ( iPollack et al.lll996l ; iKenyon fc Bromleyl 120091 ). How- 



ever, identifying exactly how dust and gas dissipate is complicated due to the difficulty in 
identifying objects clearly in the transition of losing their surrounding material. 
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One of the mos t common methods, dating back to IRAS (e.g. iStrom et all 11989 
Skrutskie et al.lll990l ). of identifying candidate transitional systems between the classical 
and debris disk evolutionary stages is through mid-IR spectral energy distributions (SEDs). 
Dust growth, sedimentation and removal is expected to result in a deficit of infrared flux 
as material is accreted or dispersed. This deficit has become the defining characteristic of 
transitional disks. For some disks, flux deficits are seen at all wavelengths, suggesting a 
gradual dissipation of mass for all disk radii. In other flux deficit is only seen at 

short wavelengths, indicating that the outer disk remains massive and optically thick while 
the inner disk lacks small dust grains. The presence of both types of transitional disks sug- 
gests that the transition happen s via different paths from the gas-rich , optically-thick stage 



to gas-poor, optically-thin disks (ICieza et al.ll2007t ICurrie et a 



with transitional disks is thought to be 5-25% (ILada et al. 



Sicilia-Aguilar et al.l l2008t iDahm fc Carpenter! 120091 ; ICurrie et al 



2009). The fraction of stars 



2006; 



Hernandez et a 



2009 



Kim et al. 



2007 



2009|). 



The small numbers indicate that the evolutionary path through a transitional disk is either 
uncommon or rapid. However, a consensus in nomenclature is still lacking: even within a 
single cluster, t he calculated fraction of transitional disks can ra nge from 10-50 %, depending 
on definitions (jErcolano et al.ll2009i ; ISicilia-Aguilar et al.ll2008l ). 



Transition disks with a short wavelength deficit are particularly interesting because they 
are potential tracers of planet formation. This deficit arises from the absence of hot small dust 
grains close to the star resulting in flux coming solely from the stellar photosphere, rather 
than disk surface emission. This deficit may be the result of grain growth and sedimentation 
but also of complex interactions between a nascent protoplanet and the surrounding disk. 
The flux deficit from such inner holes results in a depressed SED at wavelengths less than 
~15 jum, while retaining typical fluxes at longer wavelengths from the outer disk. The 
short wavelength flux deficit arises from the absence of hot small dust grains close to the 
star resulting in flux coming solely from the stellar photosphere, rather than disk surface 
emission. We define such disks as "cold disks" to differentiate from other types of transitional 
systems (e.g. anemic or homologously depleted, defined as disks with an overall lower mass 
at all radii) and to emphasize that this is an observational characteristic not necessarily tied 
to evolution. 

SEDs, relying heavily on mid-IR spe ctroscopy, are the tool currently most widely used to 



Calvet et al. 
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Forrest et al. 
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Brown et al. 



20071 . hereafter B07). Additi onal submillimeter observations of the most massive disks have 



direc t ly imaged large hole s flPietu et al.l 120061 : Hughes et al.l 120071 . 120091 : Brown et al.l 12008 



20091 : Andrews et al.l l2009h . The hole sizes match reasonably well with estimates from SED 
modeling, suggesting that current interpretation and modeling of the SEDs is correct. Within 
the category of cold disks, significant differences are seen in hole sizes (from 1-50 AU), the 
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presenc e or absence of ga s and dust within th e hole, and the presence or absence of ac- 
cretion dNaiita et alJ lioOtt lEspaillat et"alll2007l : IPontoppidan et al-lbood : ISalvk et all 12009 



Kim et al.ll2009l ). Detailed comparisons of these properties to the predictions of each of the 
scenarios listed above can help to discrimi nate which process es control the observed disk 



evolution. However, statistical studies (e.g. iNajita et al.l 120071 ) are limited by small sample 



sizes. Spitzer mapping efforts should include a significant number of unidentified cold disks 
and other types of transitional disks. However, accurately identifying these cold disks out of 
the large numbers of stars present in the maps is difficult with just broad band photometry. 
Reliable criteria for identifying disks with inner holes from Spitzer photometry, which lack 
any information on fluxes in the crucial region between 8 and 24 /im, are needed to generate 
a large sample for statistical purposes. 

This paper pre sents a new large sample of cold disks identified from 'Cores to Disks' data 
(jEvans et al.ll2009l ) and analyzes simultaneously follow-up observations of the sample with 
the IRS and MIPS instruments together with optical spectroscopy and millimeter continuum 
observations. It is organized as follows: § [2] describes the target selection, observations and 
data reduction of the different data sets, including the Spitzer, optical spectroscopy and 
millimeter observations. The results are presented in § [3] and § HJ describing, in order, the 

[31 disk masses in 



stellar parameters in § [31 disk masses in § I4.1[ the disk parameters in § 14.21 and the dust 
mineralogy in § 14. 31 for the whole sample. Notes on individual sources are given in § 14.41 The 
combined analysis is presented at § El with first the presentation of a new selection criteria 
for transitional disks in § 15.11 then a description of the observational properties of our cold 
disks in § 15.2^ and finally a discussion on the possible origins of the inner holes in § 15.31 The 
conclusions of this work are given in § [6j 



Observations and data reduction 



2.1. Target selection 



The c2d Spitzer Legacy program completed a full IRAC (3.6 to 8 /im) and MIPS (24 
to 160 /im) survey of five nearby star-forming regions (Perseus, Chameleon II, Lupus, Ophi- 
uchu s , and Serpens) ( Evans et al. 2003 ; J0rgensen et al. 2006 ; Alcala et al. 20081 ; Merin et al. 



2008 



Harvev et alil2007h . From these maps, 1024 Young Stellar Object (YSO) candidates 



were identified usi ng the selecti o n tec hniqu es, including rem oving: extragalactic contami- 
nants, described in lHarvey et al.l (120071 ) (see lEvans et al.l 120091 for a complete description of 
this sample and a general study of cloud to cloud differences). This sample should include 
a significant number of cold disks but they must first be separated from the bulk of the 
population. 
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Cold disks were selected for IRS follow-up using the following method. Spectral types 
were taken from the literature where possible or otherwise a K7 photosphere was used to 
get an initial estimate of the SED. Sources with photospheric fluxes in at least IRAC 1 (3.6 
/im) and IRAC 2 (4.5 /im) and 24 /im excesses equal or greater than IRAC 4 (8.0 /im) fluxes 
in AFa space were selected. Sources with rising or flat SEDs from 8 fim to 24 /im were 
preferentially included. A flux lower limit of 15 mJy at 8 //m was imposed to further reduce 
extragalactic contamination and keep integration times reasonable. The initially selected 
sample was then cut slightly to reject likely contaminants and to avoid overlap with other 
IRS programs. A total of 33 sources were selected in this fashion. To further expand our 
sample, two additional candidate cold dis ks within the c2d c l ouds were includ ed from a 
survey of Weak Line T Tauri stars (WTTs, iPadgett et al.ll2006t ICieza et al.ll2007l ) based on 
their IRS spectra suggestive of inner holes. In total, 35 objects were included in our sample 
and followed-up with IRS and MIPS (§ 12. 2p . optical spectroscopy (§ 12. 3p . and millimeter 
continuum observations (§ 12.41). 



2.2. Spitzer IRS and MIPS observations 

Spectra for the 35 objects in Table [1] were obtained using the Infrared Spectrograph 
(IRS) aboard the Spitzer Space Telescope under a variety of programs. Most were observed 
as part of program 30843 (P.I.: B. Merin), which also included deep 70 /im MIPS staring 
observations of the same sample. Objects 17 and 19 were obser ved with IRS as part of the 



c2d second look observations of WTTs from lPadgett et al.l ( 120061 ) and object 25 was observed 



during the IRS Guaranteed Time Observation (GTO) although it was selected as a candidate 
from the c2d photometry. All the candidate cold disks in Serpens were observed as part of a 
complete flux-limited IRS survey in Serpens, program 30223 (P. I.: K. Pontopp i dan) and will 



also be discussed along with the rest of the Serpens sample in lOliveira et al.l (120101 ) . Table 
[1] gives the AORs and observing dates for both the IRS and MIPS-70 observations of the 
sample. 

All the IRS spectra were obtained with a combination of Short-Low (SL) and Long- 
Low (LL) modules, which provide a resolving power of R = AA/A ~ 100 and a combined 
wavelength coverage from 5.3 to 35.0 /zm. Integration times were estimated from the IRAC 
flux at 8 /im and MIPS flux at 24 /im. The observations typically ranged from 5 to 20 
minutes, resulting in signal to noise values of 30 to 150, depending on the background level 
and the brightness of the objects. Given the rising SED of the targets at long wavelengths, 
longer times were usually needed in the SL than LL observations. The observations were 
scheduled in cluster mode, grouping nearby sources with similar integration times to increase 
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the observing efficiency. Data reduction started from the Basic Calibrated Data (BCD) 
images, pipeline version S12.4.0. The processing includes bad-p ixel correction, extraction 



defringing, and order matching using the c2d analysis pipeline ( IKessler-Silacci et al.l 12006 



e also the c2d Spectroscopic Explanatory Supplement, from the SSC webpagW). ThelS 
spectra can be found in Figures [1] to [5j 

The same clustering strategy was used for the MIPS observations at 70 fim. The inte- 
gration times were set to 300 seconds, with achieved sensi tivities given in Tab le |2J The small 
maps were reduced following the procedure described in ICieza et al.l (j2008a ) and using the 
SSC median-filtered BCDs from the SSC pipeline version S16.1.0. Details of the method can 
be found in the aforementioned paper, but in short, the photometry is determined with an 
aperture of 16" radius and a sky annulus with inner and outer radii of 48" and 80", respec- 
tively. Visual inspection and correction factors suggested by the SSC were used in all cases. 
The resulting fluxes are given in Table |2j together with the other IRAC and MIPS fluxes 
from c2d. 



2.3. Optical Spectroscopy 

Optical spectra were obtained fo r 21 objects in our sa mple using the Wide Field Fi- 
bre Optical Spectrograph (WYFFOS, iBingham et al.lll994l ) on the 4.2m William Herschel 
Telescope (WHT), the Intermediate dispersion Spectrograph and Imaging System (ISIS) at 
the WHT, the Intermediate Dispersion Spectrograph (IDS) at the Isaac Newton Telescope 
(INT, see Carter et al. 19940 for a description of last three instruments), the Calar Alto 
Faint Object Spectrograph (CAFOS) at the 2.2m Calar Alto Telescope, the Double Spec- 
trograph (DBSP, Oke & Gunn 1982) at the 5m Hale Telescope at Palomar Observatory, and 
with the R-C Spectrograph on the 1.5m CTIO telescope. These observations typically cover 
~ 5500 — 8500 A with a resolving power of ~ 1000 — 5000. Table [3] lists the instrument 
setups and observation log. 

Data were reduced and spectra were extracted using standard methods in IDL and 
IRAF. No telluric correction or flux calibration was performed. For WYFFOS, a fiber-fed 
multi-object spectrogr aph, sky subtractio n was obtained from fibers that were placed on 
blank sky regions (see loiiveira et aL 2009 for details). 



http : //ssc . spitzer . caltech. edu/legacy/c2dhi story .html 



http : //www. ing. iac . es/Astronomy/observing/manuals/man_wht .html 
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2.4. IRAM-30 millimeter observations 

A subsample of 17 of the northern targets was observed with the I RAM-30m telescope 



in Pico Veleta. The observations were carried out with MAMBO-2 ( IKreysa et al.l 120021 ) 
mounted on the IRAM-30m Telescope during the 2007 winter and summer pool sessions. 
MAMBO-2 is a 117-bolometer array with a half-power spectral bandwidth of 80 GHz centered 
on ~ 250 GHz (1. 2 mm), yiel ding a beam size of 11". The data were analysed with the 



MOPSI software (jZylkal Il998l ). The flux calibration was performed by observing either 
Mars or Uranus to determine the flux conversion factor. For each channel the sky noise 
was subtracted by computing the weighted mean of the signals from the surrounding six 
channels. 

All Perseus sources from our sample (see Tabled]) were observed, except numbers 3, 6, 
8 and 11. We also observed sources 20, 21 and 22 in Ophiuchus and sources 33, 34 and 35 in 
Serpens. The ON/OFF observing mode was used with a throw of 35" and each source was 
observed until a > 3a detection was obtained, or until an rms of 0.7 mJy (whenever possible) 
was reached. The integration times ranged from 10 to 40 minutes. In general, the sky noise 
of the array was low (< 80 mJy), and was never larger than ~ 100 mJy. Ten sources were 
detected with fluxes between 1 and 10 mJy, while the other 7 have upper limits between 1 
and 7 mJy. These fluxes can be found in column 8 of Table [21 



3. Stellar Properties 



Stellar properties must be determined to understand and interpret data on the sur- 
rounding disks. Optical spectra are used here to measure spectral types and to determine 
whether accretion is ongoing. We obtained optical spectra for 21 of the 35 objects in our 
sample. Of the other 14 stars, 11 have published spectral types and accretion properties. For 
the remaining three sources, spectral types are estimated from an SED analysis of broad- 
band photometry and accretion properties are left uncertain. The SED fits are also used to 
measure the extinction to each source. Results are presented in Table HJ 

For the 21 objects for which we obtained optical spectra, spectral types of G and K stars 
were assigned by finding a best match to the depth of photosp heric absorption features in 



template spectra obtained from the EXPORT spectral library (IMora et al. 



20011 ). Spectral 



types of M-dwarfs were assigned by fin ding the best match to the TiO absorption bands from 
template spectra (IMontes et al.lll997l ). The uncertainty in SpT is typically ~ 3 subclasses 
for G and K stars and 1 subclass for M stars. Object 6 is classified as a continuum object 
because no spectral features were detected in the low-resolution spectrum. Any possible 
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optical veiling was not considered for the other objects in the sample. A more detailed 
description o f our spectral typing and the stellar properties of all the Serpens objects are 
presented by lOliveira et al.l ( 120091 ). Literature values for SpT are adopted for several objects 
that we did not observe and for several objects where the literature SpT is more reliable 
than the SpT from our spe ctra. Spectral types are conver ted to T P jf using the scales in 
Kenyon fc Hartmannl (119951 ) for objects earlier than MO and lLuhman et al.l (120031 1 for those 
later or equal to MO, to account for the lower surface gravity atmosphers of mid M stars. 

The extinction and luminosity for each object were calculated using SED fits to broad- 
band optical and near-IR photometry (see Table [5] for photometry and § 14.21 for a discussion 
on the assumed extinction law). Assumed distances are listed at the bottom of Table HI 
NEXTGEN models of stellar photospheres (IHauschildt et al.lll999l ) are used as templates 
for the broadband emission. The 2MASS and any optical photometry are the primary 
constraints on extinction and stellar luminosity. For the 3 objects that lack optical spectra 
(identified with a "2" in Table H]), the SED fits are also used to estimate SpT. However, these 
SpT estimates have large uncertainties because they are based on JHK colors, which are not 
very sensitive to spectral types from mid-K through M (lLeggettlll992l ; iKenyon fc Hartmann 
1995h . 



Stellar ages and masses are estimated b y comparing t h e stel lar luminosity and effective 
temperature to pre-main sequ ence tracks of iBaraffe et al.l (119981 ) for low-mass stars (< 1.4 
Mq, mixing length = 1) and ISiess et al.l (120001 ) for higher-mass stars. Since uncertainties 
in stellar age are large they are not tabulated here and we focus the analysis on the better 
constrained stellar masses. 



3.1. Edge-on Disks 



In the process of determining the stellar properties, it became clear that a fraction 
of the sample sources were actually edge-on disks. Inferring disk structure from an SED 
analysis is difficult when the disk is viewed edge-on and occults the star. In these cases, 
optic al and near-IR light from the star and inner disk are only seen in reflected emission 



[e.g. iPadgett et al.lll999l ). The strength of such emission depends on the precise viewing 
angle and the disk flaring, but is always much fainter than it would be if the disk were 
viewed without the large absorption column from the edge-on disk. Emission from outflows, 
seen prominently in optical f orbidden lines, are not occu lted by the disk and can therefore 
have large equivalent widths (jWhite fc Hillenbrand! 120041 ) . 



Based on the existing data on our sources, we establish four criteria to determine whether 
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a disk is viewed edge-on: (1) a low ratio of photospheric luminosity to main-sequence lumi- 
nosity, assessed at 10 8 Myr from Baraffe et al. (1998); (2) large equivalent widths in optical 
forbidden lines; (3) strong silicate and/or ice absorption in the IRS spectra; and (4) high 
infrared luminosity compared to stellar luminosity. Any one of these criteria can be degen- 
erate between edge-on disks and other explanations, due to potential errors in distance and 
extinctions and the absence of accretion/outflow activity. However, the presence of at least 
two of these four criteria should establish in most cases that the disk is edge-on. Eight of 
the sources were selected as edge-on disks and the reasons are marked in Table O No further 
analysis was done on these disks as the stellar properties are too uncertain. 



3.2. Accretion Properties 



The presence or absence of accretion can be assessed from the strength and shape 
of emission in several optical lines, most prominently Ha. Though some Ha emission is 
produced in chromospheres of late-type stars, accretion typically g enerates Ha emi ssion 
with much larger equivalent widths and broader spectral profiles (see lHartmannll200ll for a 
comprehensive description of accretion phenomenology and diagnostics). A 10 A cutoff in 
equivalent width has histo rically been used to identify accretion, though this rough cutoff 
depends on spectral type (IWhite fc Basril 120031 ; iFang et al.l l2009l ). More recently, the Ha 
10% width (defined as the full width of the Ha line a t 10% of the peak flux) has been used 



to discriminate between accretors and non-accretors (IWhite fc Basril 120031 ; iMuzerolle et al. 



20031 ; IJayawardhana et al.ll2006l ). Accretion/outflow a ctivity can also produce asymmetries 



and absorption components within the line profile (e.g., Reipurth et al. 19961. 



Muzerol 



20031 . Kurosawa et al.l 120061 ) and variable line fluxes and shapes (e.g., iNguyen et al 



e et al. 



20091 ). 



For any such criteria, a small accretion rate relative to the stellar luminosity is not detectable. 

To classify the accretion properties of stars in our optical sample, a Gaussian profile is 
fit to the Ha line. We then measure the equivalent width and, when possible, the Ha 10% 
width from the Gaussian profile, after accou nting for the resoluti on of each observation (see 
a more complete discussion of this method in bliveira et aJbopgh . Uncertainties in Ha 10% 
widths range from 50-150 kms -1 . The Ha equivalent width and 10% width are then used 
together with the spectr al type to assess wheth er ac cretion is ongoing , adhering to the loose 
guidelines suggested by IWhite &: Basril (120031 ) and IFang et al.l (120091 ) . Typically, a star is 
classified an accretor if the Ha 10% width is > 300 kms" 1 and the Ha equivalent width is 
> 3 A for an early K star, > 10 A for late K star and > 20 A for an M star. 



We classify 21 stars in our sample as accretors and 8 stars as non-accretors, with 6 cases 
lacking information. In many cases an object is clearly accreting. As explained below, some 
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classifications are ambiguous both because some accretors and non-accretors have similar Ha 
10% and equivalent widths and because the uncertainties in Ha 10% width are large. Objects 
23 and 24 have weak Ha equivalent widths but broad Ha line shapes. Such observations 
could be explained if accretion is present with strong abso rption, either from o utflows or the 
accretion flow itself, that suppresses the Ha emission (e.g. iReipurth et al.lll996l ). if the star is 
an unresolved binary with only one accreting compone nt, or if the stars are chromospherically 
active with large v sin % (jJayawardhana et al.ll2006l ). The highest-resolution spectrum of 
object 24 shows absorption in Ha, while two of the three show [O I] in emission. This star 
is therefore classified as an accretor. On the other hand, object 23 shows no other evidence 
of accretion and is classified as a non-accretor. Object 18 is left unclassified because the 
Ha equivalent width is ambiguous between accretion and non-accretion. The Ha emission 
from DoAr 2 1 (object 25) has y aried on long timescales, but with no definitive evidence 
for accretion (j Jensen et al.l 120091 ) . Temporal variations in Ha also introduce uncertainty in 



classifying the presence or absence of accretion. A mass- dependent bias likely exists in this 
classification scheme because the stronger photospheric continuum emission and spectral 
type degeneracy for stars with early spectral types can mask Ha emission produced by 
weak-to-moderate accretion onto the star. Based on these uncertainties, we suspect a small 
number of objects in our sample (estimated at <10%) might be misclassified as accretors or 
non-accretors. 

To quantify accr etion, we convert t he Ha 10% widths to accretion rate using the re- 
lationship derived by iNatta et al.l (120041) . Methodolo gical differences in measuring the Ha 
10% width between this work and INatta et al.l (120041 ) likely leads to an order-of-magnitude 
uncertainty in the listed accretion rates. For late-K and early M-dwarfs, which comprise the 
bulk of our sample, the Ha 10% width is not sensitive to accretion rates < 10~ 10 M yr" 1 . 
For earlier spectral types, the sensitivity worsens to larger accretion rates. Table [4] gives the 
Ha equivalent widths, 10 % widths and mass accretion rates derived from this work and 
Table M gives the final classification of each source. 



4. Results: Disk and Disk Hole Parameters 



4.1. Disk mass 



Millimeter continuum fluxes can be used to calculate disk masses assuming that the 
millimeter emission is optically thin and a gas to dust mass ratio of 100. The disk masses 
can then be used in turn t o con strain SED models. The disk masses, MD is k, are calculated 
following iBeckwith et al.l (Il990l ) and with standard assumptions and parameters. In the 
Rayleigh- Jeans limit, 
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Mi 



Disk 



7 # C 2 

V 2A; < T>v 2 k, 



where is the flux at frequency v, here 230 GHz, d is the distance to the star (see 
Table H]), k is the Boltzmann constant, and k v is the opacity. It is assumed that the emission 
arises primarily from an approximately isothermal region with temperature, <T>. Average 
disk temperatures in the outer disk range from 10 to 50 K and here we have u sed 30 K. One 
of the major uncertainties is the opacity, k v . We follow iBeckwith et al.l (119901 ) and adopt k 
= 0.02 cm 2 g _1 anchored at u = 230 GHz, the frequency of our observations. 

Disk masses are listed in Table [7] for the ten detected sources and the seven with upper 
limits. The majority (80%) of the disks in our sample have disk masses of less than 2 x 
10 -3 M and the sample has an average Moisk of 4.5l 1 ; 4 x 10~ 3 M . These are relatively 
low ma ss disks compared to many millimeter studies, such as the large sample of T Tauri 
disks in lAndrews fc Williams! (120071 ) which have an average mass of 10 -2 M using the same 
assumptions. 



4.2. Disk hole parameters from SED modeling 



In order to determine the presence or absence of a hole, the full SED is needed to 
determine the disk structure. The SEDs for the stars in the sample are organized by type 
with edge-on disks (Figure [6]), disks with holes (Figure [7] and [8]), and disks without holes 
(Figure |9]). The solid dots are dereddened fluxes, which contains the 2MASS JHK fluxes, 
the 4 IRAC fluxes from 3.6 to 8 /xm and the MIPS fluxes at 24 and 70 /im for all sources, 
plus optical or 1.3mm fluxes for a smaller subsample of the objects. In some cases, 70 /xm 
photometry produced only upper limits, often accompanied by non-detections at 1.3mm, 
and thus severely limits information about the outer disk. The solid dots are the dereddened 
fluxes, the solid line is the IRS spectrum binned at a resolution of 1 /xm, the grey dotted line 
is the stellar Kurucz model of the central star. The differences between the Kurucz and the 
NEXTGEN stellar models used above are negligible for the purposes of disk characterization. 

The data were dereddene d using the Ay values list ed in Tableland t he ex tinction law by 
Weingartner fc Draine ( 2001 ) with Ry=5.5. Recently, Evans et alT ( 2009 ) and Chapman et al. 
( 120091 ) have suggested that this larger value of Ry is more suitable for star forming clouds 
than the typical value of Ry=3.1, which is based on models and observations of the diffuse 
interstellar medium. The higher Ry value reflects larger maximum dust grain sizes and 
results in stronger dereddening effects throughout the mid-IR. Dereddened fluxes are 2-4 
times higher at 24 /xm and Ay >10 than with Ry=3.1. However, observational measures 
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of extinction indicate that the large amount of reddening is in agree ment with observations 



and may even under estimate the actual effect in star forming regions (jChapman et al.l 12009 
Flaherty et aDl2007h . 



As the extinction law is not flat throughout the mid-IR, errors in the mid-IR extinction 
law could potentially artificially create or remove SED deficits mimicking inner holes. Cor- 
rections are much smaller at low Ay, so effects only become significant for Ay > 10. Based 
on observations, the extinction law used here is conservative with some cold disks being po- 
tentially missed due to inadequate dereddening around 20 /xm. With such large grain sizes 
in the absorbing dust, silicate grain opacities are expected to have strong effects on the 10 
fim silicate feature. In a few dereddened SEDs, strong silicate features are seen which are 
not present in the uncorrected spectra. The effects of extinction on the silicate features are 
complex and dependent on extinction and line-of- sight ice and du st composition, which may 
not be uniform throughout star forming regions (jMcClurell2009l ). Due to this uncertainty, 
all work on the silicate feature in this paper has been done on spectra which have not been 
de reddened, with the c aveat that sources with high Ay have large errors (see also discussion 



in 



Oliveira et al.ll2010h . 



The 2-D radiative transfer model RADMC (IDullemond fc Dominikl|2004l ). as modified to 
include a density reduction simulating a hole, was used to model the SEDs. The model has a 
large number of input parameters and we have fixed as many as possible. The model assumes 
a passive disk, which merely reprocesses the stellar radiation field and does not account for 
accretion luminosity and heating. This is justified since the disk to stellar luminosity ratios 
of the majority of the disks in this study are below 0.4 (see Table [7]). Further discussion on 
the disk model dependency of the results can be found in § 15.21 Stellar input parameters 
are mass, temperature and luminosity which were determined from optical spectroscopy and 
photospheric SED fits (see Table H]). A stellar Kurucz model is then used for the stellar 
radiation field. In terms of disk structure, the disk is assumed to be flared such that the 



surface height, H, varies with radius, R, as H/R oc R 2 ^ 7 , as in IChiang fc Goldreichl ( 19971 ). 
The models do not calculate the hydrostatic equilibrium self-consistently allowing the dust 
distribution to differ from the gas distribution (expected to be in hydrostatic equilibrium) 
due to effects such as dust settling. The pressure scale height is anchored at the outer disk 
edge, in this case 200 AU. In some cases, some degree of settling appears to have taken 
place. We present this number as a fraction of the scale height predicted from hydrostatic 
equilibrium, assuming a midplane temperature of 20 K. I n cases where the dis k has settled, 
a flaring index of 2/7 (0.29) is likely too large (see e.g. Andrews et al.ll2009l where flaring 
indices range from 0.04 to 0.26). Disk model SEDs with large flaring indices are similar to 
cold disk SEDs in having more flux at longer wavelengths. By using a large flaring index as 
our default, we are more conservative when reporting the presence of a hole. Surface density 
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is a power law of radius with an index of -1 inside 200 AU and then steepens to an index of 
-12 outside the 200 AU outer disk edge. The dust composition is set to a silicatexarbon ratio 
of 3:1 with only amorphous, as o ppose d to crystalline, silicate included. Silicate opacities 



were taken from iBeckwith et al.l (119901 ) which provides opacities into the millimeter. Disk 
masses were taken from 1.3mm observations where possible (see Table [7j) and left as a free 
parameter otherwise. Disk mass and disk settling are highly degenerate parameters in the 
absence of (sub)millimeter fluxes, so when the disk mass is unknown, the scale height is left 
at the hydrostatic value (see note on exception for sources 18, 19 and 24 in Section f4.4j) . 
As a result, model-derived disk masses tend to be systematically small and together with 
observational upper limits provide a range of likely disk masses. The inner edge of the disk, 
when no hole is present, was set at approximately the radius where the dust sublimation 
temperature, 1500 K, is reached. Some cold disks in the li terature have near-infrared excesses 



most easily modelled as a gap rather than a hole (e.g. iBrown et al.l 120071 ; lEspaillat et al. 



20081 ). However, none of the sources in this sample required warm dust within the hole. 
Those sources which did have some near-infrared excess turned out to be compatible with 
disk models without holes. The hole is represented in the model by a hole radius, i?Hoie and 
a density reduction, set here to a factor of 10~ 6 , providing an almost entirely clean hole. 
Thus, the four free parameters are -RhoIc, degree of settling, disk inclination and disk mass 
when there are no measurements. 

A two step procedure was used to find the best model fits. First, we ran a large grid 
of disk models with a broad range of disk and hole parameters and a limited number of 
stellar templates. Then, we ran a finer grid of disk and hole parameters, determined from 
the larger grid, for each individual star using the stellar properties in Table HI A simple 
X 2 minimization was performed at each step between the dereddened SEDs and the model 
SEDs and the resulting best fit is shown as a red dash-dot line in the SED plots (Figures [7] 
and[U]). Errors on the hole radii are given in Table [7] based on fits with up to 20% deviation 
in x 2 - 

Fifteen of the 25 modelled stars require inner holes while the remaining 10 require no 
hole to fit the SED. The 10 sources not modelled are the 8 edge-on disks and 2 sources, 25 
and 27, which have strong extended 70 /xm emission and cannot be fit with disk models (see 
section |4~41 for further details). Hole radii range from 1 AU to 55 AU. Holes smaller than 
1 AU cannot be firmly identified. The range of hole radii can be seen in histogram form in 
Figure [TTJ 

In order to confirm our results above and identify the less model-depen dent part of our 



result s, we also fitted our sample SEDs with the on-line SED fitting tool by iRobitaille et al. 



(120071 ). This system offers the possibility of fitting young star SEDs with a precomputed grid 
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of 200,000 synthetic SEDs from stars, disks and envelopes with a broad range of physical 
parameters. In order to make these fits, we observed photometry from Table [2] and |5] along 
with a binned IRS spectrum in from 6 to 33 /im with steps of 1 /im. An aperture of 4.5 arcsec, 
a distance range from 0.1 to 0.33 kpc and a Ay range from to 30 mags were used in all fits. 
For the objects for which the spectral types are known, we select the best fitting model with 
an effective temperature less than 500 K away from the target's one. The final best-fit models 
used are good representation of the median properties of the best fitting models with the 
appropriate effective temperatures but should not be considered as definitively constrained 
parameters. 

Figures [TJ and |8] show the Robitaille et al. best-fit models for the objects in our sample 
as thin black dashed lines and Table [7] shows their corresponding model ID's and the inner 
disk radii in the cases where it is clearly distinguishable from the dust destruction radius 
disk inner radius. In order to overplot the Robitaille et al. models, we corrected them with 
the extinction law given in the web-fitting tool and using the total Ay value given by the 
best-fit model and normalized it to the i^MAss-band dereddened flux. The Ay values of 
the best-fit models used are in general in agreement with the Ay values obtained from the 
spectral type and optical photometry. However, small discrepancies plus small differences 
in the standard extinction laws used to deredden our photometry and the Robitaille models 
can explain the discrepancies between the Robitaille models and optical fluxes in some cases. 

All but source 18 of the cold disks in Figure [7J also have inner holes with the Robitaille 
et al. models. The only exception, source 18, is a border-line case with an extremely small 
inner cavity of 0.03 AU in the Robitaille et al. best-fit model. On the other hand, almost 
all the objects in Figure [9j which have no hole according to the RADMC models, show 
inner cavities smaller than 2 AU in the Robitaille models. The only exception is source 15, 
which is also among the flattest disk in the group. This indicates that in general terms, the 
classification of objects as disks with or without inner holes is model independent as long as 
the disk hole is larger than 2 AU. 

Figure ITUl compa res the disk inner holes estimated both with the RADMC and the 



Robitaille et all ( 120071 ) disk model grids. The plot does not include the only outlier (Source 
18). For the rest of the sample, the values from both sets of models are correlated, but with 
a very large scatter. It must be noted that this comparison is not straight forward as both 
sets of models have different treatments of the dust opacities, disk flaring angle, and inner 
hole structure. Also, the Robitaille et al. on-line tool does not allow the use of a priori 
information like the lack of an envelope, or the effective temperature of the central star, 
both of which have obviously great impact on the final result of the fitting. Pursuing a more 
detailed comparison between both sets of disk models is beyond the scope of this paper but 
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this comparison shows that the derived disk hole radius is model independent to within a 
factor of 2-3 generally. We find the RADMC models more reliable due to the greater control 
over the input parameters and the resulting better fits. The RADMC disk radii are therefore 
used in all further comparisons. 



4.3. Dust composition 



A wide variety of mid-infrared spectral features have been det ected in spectra of disks 
aroun d low- and interm ediate mass young stars u s ing ISO (e.g 



2004) and Spitzer (e . g.. iKessler-Si 



acci et al 



2006; 



I Eur 



■ 


Acke & van den Ancker 


2006 




Lahuis et al. 


2007; 



Bouwman et al.ll2008l : IWatson et al.ll2009l : loiofsson et al.ll2009l ). These features c 
diagnostics of physical processes in disks such as grain growth, fragmentation, crystallization, 
flaring, and UV or X-ray illumination. An interesting question is whether the spectral 
features of cold disks differ significantly from those of normal T Tauri disks, and if so, what 
that tells us about the disk structure and evolution. 

Table [8] identifies the presence or absence of the 10 and 20 /im silicate emission features, 
PAHs and crystalline silicate features, where "Y" means a detection, "N" a non-detection 
and "T" a tentative detection. The spectra presented in Figures [TJ to \5\ clearly show a large 
variety of features in our sample of cold disks. 



4-3.1. Silicates 



The silicate 10 and 20 /im features are detected in most cases; no ne of the sources 



show the feature l ess spectra commonly seen for older (debris) disks (e.g., I Jura et al.ll2004l ; 
Chen et al.l 120061 ; ICarpenter et al.ll2009l ). However, only a few sources have pristine ISM- 
like 10 /im silicate features, characterized by a strongly peaked profile centered at 9.8 /im. 
This is illustrated in Figure [T3| in which the S11.3/S9.8 fim flu x ratio is plotted versus th e 
strength of the feature above the continuum S^™, as defined in lKessler-Silacci et al.l ( 20061 ). 
The different quantities 5*11.3, £9.8 and S^£™ are calculated after normalizing the observed 
feature as follows S u = 1 + (F u — F UjCOnt )/(< F U)Cont >), where F ucont is the local continuum 
and < -F^cont > the mean value of this continuum. The local continuum is estimated using 
a second degree polynomial, normalized at 6.8-7.5 /mi, 12.5-13.5 /im and 30-36 /im. The 
fluxes are measured as an averaged flux of the continuum-subtracted spectra in a ±0.1 /im 
interval at each wavelength. 



As shown by Ivan Boekel et al.l (120051 ) and lOlofsson et al.l (120091 ) , sources in which the 
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grain size distribution is dominated by large /xm-sized, or sources with size distributions 
much flatter than the MRN size distribution, grains appear in the upper left corner of 
the plot whereas those dominated by smaller ISM-type grains are located in the lower- 
right part. For compariso n, th e silicate profile para meters of the c2d sample studied by 
Kessler-Silacci et al.l (120061 ) and lOlofsson et al.l ( 120091 ) are overplotted (grey dots). It is seen 
that the cold disks fall mostly in the upper left part of the figure, indicating that the dust 
in these disks has already undergone non-negligible grain growth. No significant difference 
was found in the 10 /xm shape for accreting or non-accreting objects. This conclusion is not 
affected by crystallization or sedimentation of dust toward the midplane, which causes some 
spread in the relation but can not reproduce the observe d trend, as discussed extensively in 
Dullemond fc Dominikl (hoosh and blofsson et all (l2009h . 



Some confirmed disks with inner holes show crystalline silicate features, either the 



forsterite feature at 33 /im and/or the 23 and 28 /xm complexes (see lOlofsson et al.l 12009 



for definition of these features). The fraction is between 33% and 60% (5-9/15) with the 
range due to some tentative detections limited by the low signal-to-noise ratio in some spec- 
tra. A prominent example of a highly crystalline spectrum is that of object 10, SSTc2d 
J034227. 1+314433, a cold disk in Perseus (see Figure [5]). This detection frequency of crys- 
talline material is comparable to the 55% f ound for the c 2 d sam ple of normal CTTS disks 
based on long-wavelength features, whereas I Wat son et al.l (120091 ) claim an even higher frac- 
tion of 94% for normal disks in Taurus (see discussion in Olofsson et al. 2009). 



4.3.2. PAHs 



The detection of PAHs is based on at least two of the features at 6.2, 7.7, 8.6, 11.2, 12.8 



and 16.4 /im (IKoike et al.lll993l : iGeers et al.ll2006l : iTielend 120081 ). Overall, PAH features are 
detected in 13% (2/15) of the cold disks in our sample. This is similar to th e fraction of 



11-14% found by lGeers et al.l (120061 ) for a sample of normal CTTS disks (see also lGeers et al. 
2007bl ). This low detection fraction not only reflects the fact that late- type stars have less 



UV and optical photons to excite the PAHs but also implies a PAH abundance in disks 
that is typically a factor of 10-100 lower than in the normal interstellar medium, taken to 
be 5 x 10 ~ 7 with respect to total hydrogen for PAH species with 100 carbon atoms (e.g., 
C100H24) (IGeers et al.l l2006l ). Both sources 25 and 27, which have long wavelength fluxes 
composed of extended material, show exceptionally strong PAH features. Both have low 
continuum at mid-infrared wavelengths increasing the line-to-continuum contrast. In the 
case of source 27, even the 16.4 /im band is clearly revealed (see Fig. IT2l) . 



The strong features seen in cold disks do not necessarily imply higher PAH abundances 
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than in normal CTTS disks, however. iDullemond et all (120071 ) model the PAH emission 
in disks with varying degrees of sedimentation, as characterized by the a parameter of 
turbulence. In the case of low turbulence, a < 10~ 4 , the big grains quickly sediment to the 
midplane but the PAHs stay in the upper layers, boosting the feature-to-continuum ratio by 
a factor of 2-10. For high turbulence, a ~ 10~ 2 , the small and big grains stay well mixed 
in the upper l ayers, resulting in t he typically weak or abs ent PAHs features consistent with 
observations (IGeers et all 120061 ) . IDullemond et all ( 120071 ) did not consider disks with inner 
holes, but qualitatively the effect is similar to their low turbulence case in which the big 
grains are removed and the features strengths are boosted. 

In contrast with the silicates, the PAH emission arises not only from the inner disk but 
also from the outer disk surface exposed to optical and UV radiation. Indeed, spatially ex- 



tended PAH emission ha s been detected out to ~60 AU or more for some disks (IHabart et al. 



20061 . IGeers et al.l l2007aJ ) . T wo cold disk s within the c2d clouds, SR 21 and T Cha, are in- 
cluded in the observations of IGeers et al.l ( 12007bl ) using VLT-ISAAC and VLT-VISIR to get 
higher spatial resolution compared with Spitzer. In both cases, the emission is found to be 
spatially unresolved, with limiting spatial extents (FWHM) of 19 and 13 AU, respectively. 



Thes e limits are comparable to the values of i?Hoic of 18 and 15 AU, respectively (IBrown et al. 



20071 . Table [7]), indicating that the PAHs are located primarily inside the gaps. For the case 
of SR 21, the location of the ga s inside the gap has been pinpointed to a ring at ~ 7 AU 
radius ( iPontoppidan et al.ll2008l ). proving that the gas emission is indeed coming from inside 
the gap and not from the outer edge or wall of the gap. Since PAHs are likely coupled to 
the gas, PAH emission, when present, serves as confirmation of the presence of gas inside 
the gap. Unfortunately we see no correlation between PAHs and accretion, both of which 
should trace gas within the holes. This is likely due to the dependence of PAH emission on 
strong UV flux which is lacking from late-type stars. 



4.4. Notes on individual sources 

Sources 10, 12, 14, 20, 22, 23, 29 and 30 have neither 70 fim nor 1.3 mm fluxes available. 
This makes it very difficult to distinguish between a low mass disk and a more massive disk 
with an inner hole. The difference between the two lies in the presence of a substantial outer 
disk. The 70 /im limits sometimes hint at the absence of such an outer disk but the limits 
are not always very stringent. Unless there is compelling evidence from the long wavelength 
IRS spectra, we have taken the conservative view that these are likely low mass disks and 
can be explained without needing to invoke an inner hole. 

Two sources, 25 and 27, have extremely strong, extended 70 /iin emission and do not 
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appear to be disks. In both cases, it is likely that surrounding material has contaminated 
the long wavelength photometry, producing these strange SEDs (see Figure [6]). Source 25, 
DoAr 21, has been better studied than many of the other sources in this sample. Resolved 
imaging of H 2 shows an e xtended ring structu re at 73-219 AU away from the star (Hogeheijde 



et al. 2009, submitted). iJensen et al.l (120091 ) have speculated that much of the flux in long 
wavelength unresolved photometry comes from cloud material that the star is heating. Thus, 
while the SED looks compatible with a large inner hole, this might be due to excess material 
not directly associated with the star with long wavelength fluxes being increasingly affected. 
Source 27 has a similar SED and we suspect may arise from a similar physical situation. 
Both objects show exceptionally strong PAH emission (see Figure [121 and § I4.3.2j) . Within 
a disk structure, the strong 70/im fluxes indicate large amounts of dust at 100-200 AU 
and good SED fits required unphysical input parameters. Objects with very pronounced 
70 /im excesses but small 24 fim excesses have been found to be produced by background 
contamination in the great majority of the cases (Wahhaj et al. 2009, submitted). For these 
reasons, we have removed these objects from further study. 

The 70 /im flux of source 17 is affected by strong variable background emission, clearly 
visible in the c2d MIPS-70 image. Two photometric points are given in the SED (Figure [7]) 
with the larger being the full aperture flux and the smaller a PSF aperture flux. The actual 
flux may be even smaller. The 24 /im image does not show background emission and is in 
good agreement with IRS spectra. Although there is evidence for background contamination 
in the long wavelength SED of this object, a sensible fit was possible with a disk with a hole. 

Sources 11, 12, 13, 17, and 18 have stellar temperatures lower than the Kurucz model 
grid and a blackbody is used instead. 

Sources 18, 19 and 24 form a distinct subset of our sample. These three disks all have 
little to no near-IR excess but very flat disks longward of 8 /im. While there does appear to 
be a discontinuity between the inner and outer disks, the difference is not large. Despite the 
unknown disk mass, some degree of settling in the outer disk was required to reproduce the 
SED shape. 

Finally, it is worth noticing that objects 2, 13 and 33 (15% of the cold disk sample) 
show factor of 2 differences between IRAC and IRS fluxes in the overlapping wavelength 
range which are too large to be due to calibration errors. In the particular case of object 
33, which has the largest difference, only the IRS spectrum hints at a possible inner hole, 
while the IRAC photometry could be easily explained without a hole. A careful check of 
the extraction of the IRS spectra and IRAC photometry found no evidence of extended 
emission, mis-alignment in the IRS observations, or any other instrumental reason for the 
discrepancy. Indeed, only 6 objects (4%) out of the 147 IRS observed disks in the Serpens 
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IRS survey have similar flux discrepancies between IRAC and IRS (Oliveira, I., 2010, priv. 
comm.). The IRAC observations were taken approximately 2 ye ars before the sp e ctrosc opic 
ones (see Table [T]). Similar scale variability was reported in iMuzerolle et al.l ( 120091 ) for 
another transitional disk. Such intrinsic variability might be the cause of this phenomenon 
but systematics between spectra and photometry could cause complications and further 
investigation is outside the scope of this work. 



5. Discussion 
5.1. Robust Selection Criteria 

5.1.1. Spectroscopic selection criteria for cold disks 
Determining concrete selection criteria for cold disks is necessary to identify cold disks 



Brown et al. 


2007; 


Fur lan et al. 


2009) 



samples using 30/13 /im flux ratios and slopes (e.g. 
These two wavelengths bracket the rise in flux between an inner dust hole and substantial 
outer disk for holes with radii between 1 and 100 AU around stars with a representative 
range of stellar masses. The disks in this sample provide confirmation of this classification 
with 93% (14/15) of the cold disks having F 30 /Fi 3 flux ratios between 5 and 15 (see Figure 
[T41) . The only disk with F 30 /-Fi3<5 is Source 2, which has one of the smallest holes. Edge-on 
disks generally lie above i<3o/-Fi3=15 but three overlap with the cold disk region. There is 
also not a clean separation between normal and cold disks at F3o/-Fi3=5. The 30/13 /im 
flux ratio has a crude correlation with larger hole sizes having larger flux ratios, but there 
is unfortunately no clear relation between the 30/13 /im flux ratios and the slope between 
Spitzer 8 and 24 /im photometry, particularly for cold disks. Thus 8 to 24 /im slopes and 
colors cannot be used as a proxy for estimating hole size, although they can be used to select 
cold disks, as we discuss in the next section. 



5.1.2. Photometric selection criteria 



To identify the cold disk candidates out of the YSO population, we tested several com- 
binations of color-color diagrams with different Spitzer bands to determine the cleanest 
identification of disks with inner holes (cold disks). We included not only the results of the 
detailed SED fitting from § 14.21 b ut also literature in formation about cold disks in the c2d 
clouds including two sources from iBrown et al.l (120071 ) and two sources from Andrews et al. 
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( 120091 ). The selection criteria proposed here can be applied to any Spitzer IRAC and MIPS 



YSO photometric sample of any star-forming region to identify the disks with inner holes. 

The criteria are based on colors without dereddening as the extinction to many of the 
sources in the catalog is unknown. The effects of Ay=10 are plotted as an arrow on the 
diagrams. The effects are small for extinctions less than 10. Analysis of the spectral sample 
based on deredden colors provides little difference to the selection regions. 

The cold disk candidate sample is defined with the following cuts in magnitudes, orga- 
nized in two sections for two different types of objects (see top left Figure [T5|) : Region A 
selects "clean" inner holes (i.e. disks for which there is no substantial excess in any IRAC 
band and the signature of an inner opacity hole). Region B selects disks with some excess 
disk flux in the IRAC bands. The definition of the boundaries are the following: 

Region A : 0.0 < [3.6] - [8.0] < 1.1 ; 3.2 < [8.0] - [24] < 5.3 (2) 



Region B : 1.1 < [3.6] - [8.0] < 1.8; 3.2 < [8.0] - [24] < 5.3 (3) 

Region B contains the 4 literature sources plus source 2, while the majority of the new cold 
disks (13/15) lie in Region A. One possible explanation for this division is the focus of this 
work on looking for photospheri c fluxes in the nea r -IR, while disks with large holes often 



have small near-IR excesses (e.g. iBrown et al.l 120071 . lEspaillat et al.ll2007f ). 



From the overall photometric catalog, the statistics on cold disks within the five c2d 
clouds can be determined, including sources not included in the spectroscopic follow-up for 
a variety of reasons such as low fluxes. These disks are marked with small crosses in Figure 
[T5l There are 43 cold disk candidates in Region A accounting for 4% of the total disk 
population in the c2d YSO catalog. Only one of the objects in that region for which we 
performed detailed SED modeling including the IRS spectrum was a false positive. This is 
object 25, DoAr 21, which has a cold disk-like SED but the 70 /xm flux may be strongly 
contaminated by cloud material (see § 14.41 for further details). From the spectroscopically 
studied objects in Region A, 93% are cold disks. 

Region B is more complicated. While there are well studied, confirmed cold disks 
within this region, there are many objects which are not cold disks. Region B contains 181 
YSOs in the c2d catalog, but fewer sources were selected for our study from this region. 
Of this small subsample and excluding the rising sources, for which it is not possible to 
determine the presence of an inner hole in the disk, only 42% of the studied disks turned 
out to be cold disks, so we estimate a contamination on the order of 58% in this region. 
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These percentages have considerable uncertainty given the low-number statistics. Region B 
includes the disks WSB 60 and DoAr 44, for which inner holes of 20 and 33 AU, respectively, 
have been resolved w ith the SMA but whose SEDs show no detectable signature of them 
( [Andrews et al.ll2009f ). This implies that selecting purely on broadband colors will always be 
prone to uncertainty. Individual source modeling, spectroscopic studies and non-SED based 
searches are all needed to fully understand cold disk frequencies, particularly in this region 
of the color-color parameter space. 

If we add together the 43 objects in region A with 42% of the 181 objects in region B, 
we have a total of 119 cold disks in the sample (12% of the YSO catalog of 1024 sources). 
Therefore, the fraction of 4% of cold disks quoted above should be interpreted as safe lower 
limits to the actual population of disks with inner holes in the c2d YSO catalog. The 
analysis of this photometric sample is outside the scope of this paper and will be presented 
in a separate article. Considering the IRS flux-limited sample in Serpens where all disks 
have IRS spectra, we confirm the presence of 8 cold disks, which account for 9% of the total 
( lOliveira et al.ll2010l ). This range of cold disk frequencies of 4-10% comes with some caveats. 
Transitional disks wit h very low mass outer disks are no t included in this study (as the 
ones discussed in e.g. iNajita et al.l 120071 ; ICieza et al.l 120071 ) . but are often included in other 
transitional disk frequency statistics. Such disks may represent a later stage of cold disk 
evolution or a completely different evolutionary pathway. Disks with holes smaller than 1 
AU cannot be reliably identified. Disks with large holes but with some near-IR excess can 
also be difficult to identify from photometry. 

A final consideration is the potential bias introduced by the c2d YSO selection criteria 
themselves. The c2d YSO sample is biased against objects with photospheric colors and only 
includes sources with good detections at all IRAC bands plus at MIPS 24 /jm band. This 
obviously imposes bias which is a difficult to determine against selecting transitional disks 
in regions with very high background emission, like e.g. the IC 348 or p Oph clusters. In 
any case, the effect seems not to be a dominant one since most of the cold disk candidates 
found with the color selection proposed here are clustered in a similar manner as the rest 
of the YSOs around the dense high-background clusters, while we should see a lack of such 
objects in high background areas otherwise. 



5.1.3. Comparison with transitional disk selection criteria from other papers 



Fang et al.l ( 120091 ) also recently proposed color-color diagram cuts to select cold disks. 
Their criteria form a trapezoid in K-[5.8] vs. [8.0]-[24] color-color space defined by K-[5.8] 
> 0, [8.0]-[24] > 2.5 and the last side by K-[5.8]<(0.56+([8.0]-[24])x0.15). In Figuredgb, we 
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evaluate these criteria based on our spectrally confirmed cold disks. Most of the disks in our 
sample fall within or close to this area. Of our cold disks, 13/15 are within the boundaries 
while the remaining two are only just above the upper sloped line. However, half (6/12) of 
the disks without holes also lie within this region as well as some (3/8) of the edge-on disks. 
Cold disks can be preferentially selected by cutting higher in [8.0] -[24] at ~ 3.5 instead of 
2.5. This would exclude more of the candidates that do not have inner holes while including 
all the disks with holes and also have the benefit of having a greater separation from the 
bulk sample in the larger catalog (shown in the figure as small gray crosses). The region of 
[8.0]-[24] between 2.5 and 3.5 has the potential to include transitional disks with low mass 
outer regions, but individual disk modeling is needed to confirm the nature of these sources. 
An upper limit in [8.0] -[24] is also needed to remove the edge-on disks, with [8.0]-[24]<5 
producing a reasonably clean separation. 



Muzerolle et al.l ( 120101 ) use slopes in log \F\ vs. log A between 3.6 (im. and 5.8 /im, 
«3_5, vs 8 /im to 24/xm, a?8-24 • They propose that all transitional disks lie below -1.8 in 
«3_5 with weak excess sources between 0:3-24 = — 1.5 and 0:3-24 = an d normal cold disks 
at «8-24 >0. Our sources lying within the weak excess region were all found not to need 
holes to fit the SEDs. The normal cold disk reg ion selects similar sour ces to Region A but 



sources within Region B are not selected by the iMuzerolle et al.l ( 120101 ) criteria 



Cieza et al.l ( 120101 ) adopt a wider definition of transitional disks (e.g. they include 
objects with very small 24 fim excesses), and their selection criteria are therefore less useful 
for identifying cold disks. Transitional disks are taken to be in the quadrant defined by [3.6]- 
[24]>2 and [3.6]-[4.5]<0.25. Just over half (12/19) of the cold disks lie within the selected 
quadrant, with 7 cold disks missed by this criteria (generally those with more massive outer 
disks). Almost half (5/12) of the disks without holes lie within the region and in the exact 
same position as the disks with holes. These disks may be homologously depleted transitional 
disks. A fairly large fraction of the catalog classical disk locus lies within the selected region 
and might lead to a large false detection rate of cold disks. 



5.2. Properties of the cold disks 



One of the most obvious characteristics of this cold disk sample is the relatively small 
hole sizes (see distribution in Figure [Til. Many of the IRS selected cold disks in th e literature 



Calvet et al. 2005; Brown et al. 



20071 : lEspaillat et al.ll2007l : Kim et alJl2009h have large 



(e.g. 

holes with radii greater than 10 AU. Of the 15 cold disks examined here, only three have 
hole sizes larger than 10 AU while an additional four large hole sourc es within the clouds 
are identified in the literature (IBrown et al.l 120071 : 1 Andrews et al.l 120091 ) . This may indicate 
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that cold disks with very large holes (>10 AU) are intrinsically rarer than disks with smaller 
holes. Large hole sizes are easy to see from IRS spectra but may not be as obvious from 
photometry. One complication is that disks with large inner holes often have some near-IR 
excess making the disks difficult to identify from IRAC colors. Another possible complication 
arises from the different disk models and methods used to determine the inner disk hole radii. 
In particular, the possibility cannot be ruled out that the RADMC models, which do not 
include viscous heating due to acc retion tend to pro duce smaller inner holes than other 



accreting disk models used by e.g. ICalvet et al. 



Dn t he ot her hand, the hole radii 



determined with same models and methods by iBrown et al.l (120071 ) were confirmed by the 
direct SMA observations of three disks by Brown et al. (120091 ) . Finally, the si milarity in terms 



of inf rared SEDs of the two millimeter-discovered cold disks in Ophiuchus (jAndrews et al. 



20091 ) to classical disks means millimeter surveys, as well as infrared, would be needed to 
ensure that all cold disks were found and to robustly c onfirm this tr e nd. T he d ifferences in 



the defi nition of transitional disks between the works of iNaiita et al.l (120071 ) and ICieza et al. 
( j2008at ) were also extensively discussed in ICieza et al.l (120101 ) . 



We find a statistically significant positive correlation between disk mass and hole size (see 
Figure [TB]) . The correlation has a Pearson's correlation coefficient of 0.6 and is statistically 
significant on the 99% level including upper limits on disk mass and on the 90% level with 
only detec t ed so urces. To minimize systematic differences, disk ma sses for sources from 
Kim et al.l (120091 ) have been recalculated using the 1.3 mm fluxes in Andrews fc Williams 
( 120051 ) and the flux to disk mass conversion in Section 14.11 The dependence of hole size on 
disk mass may point to a gravitational process, with more massive disks more likely to form 
more massive planets in larger orbits . The masses of transitional disks relative to normal 
disks remains unclear. INajita et al.l ( 120071 ) finds that transitional disks have larger disk 
masses for a sample of transition al disks in Taurus, although their disk classification method 
is different. ICieza et al.l (l2008bf ) find very small disk masses for a sample of non-accreting 
W eak-Lined T Tauri stars. This trend may help explain the different results of the studies 
as INajita et al.l ( 120071 ) focused mainly on cold disks from the l iterature which generally have 



large holes and are often accreting, while ICieza et al.l ( 12008bl ) were likely more sensitive to 
disks with smaller holes and no accretion. 



Kim et al.l ( 120091 ) find a correlation between hole size and stellar mass in their transi- 



tional disk sample in Chamaeleon and Taurus. Unfortunately, our sample is too similar in 
stellar mass to accurately test this trend and there i s no correlation found from our sam- 
ple alone. However, combining with the iKim et al.l (120091 ) sample (see Figure IT7|) retains 
a correlation between the size of the inner holes and the stellar masses. The disks from 



Brown et al.l (120071 ) are in agreement with the IKim et al.l ( 120091 ) trend while the disks in this 
paper generally lie below the trend line, indicating smaller hole sizes. This may be partly 
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due to the generally lower disk masses in this sample leading to smaller hole sizes based on 
the correlation in Figure [T6l 

Perhaps surprisingly, we see no trends with L^^/L^, a measure of overall disk evolution 
(see Table [7]). L-^^/L^ measures the integrated infrared excess normalized by the stellar 
luminosity and primarily traces the total grain surface area that is reprocessing stellar light. 
As disks become more tenuous, settle and disappear, the strength of the d isk luminosity 



shoul d decline. This measure is commonly used in debris disk studies (e.g. iMorales et al. 



20091 ). There is no correlation with hole size nor are there significant differences in Ld^/L^ 
between the sources with and without holes in this sample. Therefore, if I>Disk/£* can indeed 
be used as a probe for disk flaring, then this points to a distinct lack of settling in the outer 
regions of cold disks. Alternatively, this might also imply that L^^/L^ is not especially 
sensitive to the changes in the disk flaring while the outer disk keeps being optically thick. 

The majority (9/12, 75%) of the cold disks with Ha spectra are accreting. Accretion is 
a particularly interesting diagnostic in transitional disks as it indicates material inside the 
hole. Large fractions of accreting cold disks indicate that in many cases gas must be flowing 
through t he hole as any inner gas reservoir would drain quickly. This is in contrast with the 



results by iNajita et al.l (120071 ). who report lower mass accretion rates in transitional disks 
with respect to the classical T Tauri stars in Taurus, although using a completely different 
definition of transitional disk from the one used here. 

Silicate emission probes dust properties such as grain size and crystallinity. Silicate 
emission, especially the distinct 10 /im silicate stretching band, is seen from all but one of the 
cold disks. However, the presence and intensity of the 10 /im feature decreases substantially 
for inner holes radii larger than 7-10 AU (Figure fT8l) . This is partly to be expected as 
dust-depleted inner holes result in a lack of silicate grains at the correct disk temperatures 
to produce the features. A large fraction of the cold disks show crystalline features at 
wavelengths longer than 20 /im indicating processing which likely occurred in the inner 
regions of the disks where temperature and densities are high. These regions are now largely 
depleted in the cold disks but dust mixing, fragmentation of larger bodies or alternative 
formation routes have resulted in detectable amounts remaining further out in the disk. 



5.3. The possible origins of the inner holes 



The origins of inner holes are still under debate among several theories developed to 
explain inner hole s and gaps in pro t oplanetary disks , inclu ding i) EUV photoevaporation 
of the inner disk (IClarke et al.ll200ll ; Alexander et al.ll2006l ). ii) settling and coagulation of 
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dust into large particles f lTanaka et al.l l2005r iDullemond fc Dominikl 120051). or in) by the 
dynamical pe rturbation of unresolved companions, bo th stars ( lArtymowicz fc Lubowlll994l ) 
and planets (I Guillen et al.l 120041 ; IVarniere et al.l 120061 ) . It is likely tha t all are pre s ent o n 
some level but it is not c lear which are dominant in which systems ICieza et al.l ( 120101 ); 
Sicilia-Aguilar et al.l fcoioh . 



For photoevaporation, an inner hole occurs when the photoevaporation rate driven by 
the E UV ionizing flux from the central star matches the viscous accretion rate ( lAlexander et al 
20061 ). FUV photoevaporation also plays a significant role in dissipating disks but predom- 
inantly rem oves less bound gas from the outer regions and would therefore not create cold 
disk SEDs (IGorti fc Hollenbachl |2009| ). Photoevaporation is most effective when accretion 
rates are low and would result in no gas or dust close to the star within 0.1-1 Myr, de- 
pending on assumed disk properties such as visco sity and treatment of the UV ionizing flux 
( lAlexander et al.ll2006t IGorti &: Hollenbachl 120091 ). The presence of accretion in such a large 
fraction (75%) of this sample discards photoevaporation as the primary origin. Some disks 
with inner holes do show evidence for photoevaporation based on [Ne II] o bservations but the 
impli ed mass loss rates are too low to disperse the disk in under 10 Myr ( jPascucci &: Sterzik 
20091 ). The disk masses of the remaining 25% are typically low enough that photoevaporation 
could be the cause of the inner hole. 

Accelerated grain growth in the inner regions coul d produce cold disk SEDs as dust 
grains became too big to be observed (ITanaka et al.ll2005l ). Models predict that dust particles 
grow and settl e towards the dense d isk mid-plane, where they may stick together to form 
planetesimals ( jWeidenschillingl l2000h . Growth is likely preferential in the inner disk, and 
these larger bodies will gro w and can eventually accrete a large frac tion of the surrounding 
gas to become giant planets (IBryden et al.lll999l ; IWuchterl et al.ll2000l and references therein). 
In this scenario, grain growth and settling would happen throughout the disk but with faster 
timescales in the inner region. If grain growth is accompanied by settling throughout the 
disk, decreasing Loisk/-^* with increasing hole size would be expected. We do not see this 
for our sample making this unlikely to be the dominant scenario although such a correlation 
for a subsample of the sources may be masked by short inner disk settling timescales or 
objects formed by different mechanisms. Another potential complication is that coagulation 
in the inner regions would not effect L^^/L^ in the absence of settling in the outer disk. 
Holes created via grain growth are likely to have more gradual hole edges compared to holes 
created by companions. 

Binary companions are a possible explanation for cold disk SEDs. So me candidate cold 
disks , such as CoKu Tau 4, have been later determined to be circumbinary (llreland fc Kraus 
20081 ). The issue clearly requires further high-resolution imaging of these sources (e.g. 
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Marois et al. 


2008; 


Kalas et al. 


2008; 



to confirm the possible presence of stellar or substellar-mass companions. However, this is 
time consuming for large samples and may not even be possible for more distant regions 
or planets at larger radii which require long term monitoring. Preliminary results of bi- 
nary searches show that close binary companions often result in the com plete dissipation of 



the disk and tha t circumbinary "transition" disks may not be common (IKraus et al.l 12009 



Pott et al.ll2010l ). If stellar companions carve out inner holes, accreti on may either continue 
or ce a se depending on the se paration and mass ratio of the two stars (lArtymowicz fc Lubow 
19961 ; llreland fc Krausll2008l ). On the other hand, if a planet carves out an inner hole ac- 
cretion can continue at a reduced rate ai ded by mechanisms such as the magneto- rotational 
instability Jchiang fc Murrav-Clavll2007t ). 



Inner holes created by planets remains a popular and exciting explanation. Of the 
scenarios described this is the most difficult to positively confirm or deny. Planets are 
faint so direct observation is difficult while their small masses leave weaker gravitational 
signatures on the disk than stellar companions. The inner hole radii are compatible with 
the distribution of exoplanet semi-major axis, as found in the latest exoplanet database 
as of December 2009^ (Fig. [T7|) . As the exoplanets likely formed in similar disks before 
dissipation, it is likely that there are young exo-planets in some disks at these types of radii 
but whether they are responsible for the cold disk signatures discussed in this paper remains 
to be determined. 



6. Summary 

The main results of this work are summarized here: 



• Optical spectra, 2MASS and Spitzer photometry, millimeter continuum observations 
and Spitzer/IRS 5 to 35 /im spectra of a sample of 35 cold disk candidates selected 
from c2d photometry are presented and analyzed. 

• Out of 35 objects in the initial sample, SED modeling identifie s 15 as disks with i nner 
holes, which we call "cold disks", following the c2d convention (IBrown et al.ll2007l ). Of 
the remaining sources, 10 could be modelled without holes, 8 are edge-on disks and 2 
have SEDs strongly contaminated by cloud material. 



• fottp : //www . e xoplanets . eu 
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The color cuts 0.0 < [3.6] - [8.0] < 1.1 and 3.2 < [8.0] - [24] < 5.3 identify most 
cold disks from this sample (~ 80%), in particular those with the cleanest inner holes. 
Extension of the [3.6] — [8.0] color cut to 1.8 recovers some objects with small near-IR 
excesses and large holes, but contains contamination from disks without holes. Out of 
the large c2d YSO sample, between ~ 12% of the disks are estimated to be cold disks 
based on these selection criteria. 



We ev aluated the criteria of iFang et al.l (120091 ). iMuzerolle et al.1 (120101 ) and lCieza et al. 



(120101 ) and suggest improvements based on our spectroscopic study. 



The cold disks presented here have small hole sizes, generally less than 10 AU. This 
distribution is more in agreement with exoplanet radii than the large hole sizes of most 
cold disks in the literature. 

A large fraction (75%) of the cold disks are accreting, suggesting that gas is flow- 
ing through the dust depleted hole. This large fraction of accreting disks is not in 
agreement with the dominant hole origin being photoevaporation. 

Hole size correlates with disk mass with more massive disks tending to have larger 
holes. 

The sizes of the inner holes scale linearly with the stellar mass although with a large 
spread. 

The 10 /jm silicate features in the sample show substantial grain growth. The 10 /xm 
silicate emission feature strength with respect to continuum decreases drastically for 
inner holes larger than ~ 7 AU. Some (33-60%) of the cold disks show long wavelength 
crystalline features indicating that mixing from the inner regions where crystallization 
occurs to outside the inner hole region must be efficient. Only 2 sources (~ 13%) show 
PAH emission. 



Support for this work, part of the Spitzer Space Telescope Legacy Science Program, 
was provided by NASA through Contract Numbers 1224608 for c2d, 1288664 for G03 IRS 
cold disks, 1256316 and 1230780 issued by the Jet Propulsion Laboratory, and for California 
Institute of Technology under NASA contract 1407. Astrochemistry at Leiden is supported 
by a NWO Spinoza and NOVA grant, and by the European Research Training Network "The 
Origin of Planetary Systems" (PLANETS, contract number HPRN-CT-2002-00308). 
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Facilities: Spitzer (IRAC,MIPS,IRS), CAHA 2.2/CAFOS, WHT (ISIS,WYFFOS), INT 
(IDS), Palomar/DBSP , Pico Veleta/MAMBO. 
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Table 1: Cold disk candidates and Spitzer observing log 
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Table 2. IRAC, MIPS and MAMBO 1.3 mm fluxes for the cold disk candidate sample 
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z 
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i 
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Z 


OQn i m j_i nVn 1 nn 


1 
1 


79e+01±4 


47c+00 


^ A QOn 1 no 
<4.oze-|-Uz 






24 


6 


07c+01±2.74c-01 


■i 


17c+01±2.74e-01 


3 


09e+01±1.02e-01 


2 


28e+01±1.20e-01 


1 


59e+02±9 


73e-01 


5.37c+02±7.88c+01 






25 


1 


19e+03±1.94e+01 


8 


45e+02±1.49e+01 


7 


30c+02±1.00e+01 


6 


39e+02±8.50c-(-00 


1 


34e+03±3 


42e+01 


1.20e+04±1.26e+03 






26 


9 


80c+00±5.75e-01 


8 


51e+00±5.75e-01 


7 


03e+00±3.82e-01 


5 


46c+00±2.98e-01 


4 


34e+01±4 


03e+00 


<2.87e+02* 


2.40e+01±6 


00e-01 


27 


5 


25e+01±2.97c+00 


3 


67c+01±2.97c+00 


3 


12e+01±1.68e+00 


2 


85e+01±1.68e+00 


9 


74e+00±9 


20e-01 


1.17c+03±1.01e+02 






28 


2 


13e+01±1.43e+00 


2 


93e+01±1.43e+00 


2 


25e+01±1.29e+00 


2 


27c+01±1.10e+00 


1 


02e+01±3 


72e+00 


<2.87c+02* 






29 


4 


46e+00±2.73c-01 


3 


36e+00±2.73c-01 


2 


39c+00±1.47e-01 


1 


72e+00±1.03e-01 


1 


60e+01±l 


49e+00 


<2.87e+02* 






30 


1 


36e+01±8.67e-01 


1 


62e+01±8.67c-01 


1 


51e+01±8.86e-01 


2 


56e+01±1.26e+00 


6 


05e+01±5 


59c+00 


<5.00e+01 






31 


6 


42e+00±4.61c-01 


7 


33c+00±4.61e-01 


5 


99e+00±3.46e-01 


8 


05e+00±3.94e-01 


3 


25e+01±3 


Olc+00 


1.35c+02±1.57c+01 






32 


3 


05e+01±2.00e+00 


2 


41e+01±2.00c+00 


1 


70c+01±9.68e-01 


1 


28e+01±6.61c-01 


7 


42e+01±6 


86e+00 


5.76e+01±1.31e+01 






33 


2 


lle+01±1.24e+00 


1 


99e+01±1.24e+00 


1 


56e+01±8.31c-01 


1 


28e+01±6.19c-01 


8 


llc+01±7 


49e+00 


1.26e+02±1.74e+01 


2.30e+00±6 


00e-01 


34 


2 


08e+01±1.19e+00 


1 


37e+01±1.19e+00 


1 


10e+01±5.91e-01 


6 


97e+00±3.45e-01 


1 


97e+01±l 


83e+00 


1.40e+02±1.71e+01 


6.30e+00±6 


00e-01 



Table 2 — Continued 



No. 


IRAC 3.6 


IRAC 4.5 


IRAC 5.8 


IRAC 8.0 


MIPS 24.0 


MIPS 70.0 


1.3 mm 




(mjy) 


(mjy) 


(mjy) 


(mjy) 


(mjy) 


(mjy) 


(mjy) 


35 


1.04e+01±6.24e-01 


9.28e+00±6.24e-01 


8.20e+00±4.53e-01 


1.06e+01±5.43e-01 


1.82e+01±1.69e+00 


1.76c+01±6.81e+00 


2.20e+00±9.00e-01 



* Upper limits are estimates from nearby sources due to difficulties with the local background flux. 
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Table 3. Optical spectroscopy observation log 



Date 


Telescope + 


Grism 


Wavelength 


Dispersion 


colheadID 




Instrument 




Coverage 


(A/resel) 




2006 November 22 


WHT+WYFFOS 


R600R 


5500-8500 


4 


2, 3, (4), 5, 7, (10), (11 












12, 13, 14, 15 


2006 December 13 










1,(9) 


2007 July 26 


WHT+ISIS 


R600R 


6000-8000 




24 


2008 July 22 










2 


2007 May 06 


INT+IDS 


R600R 


4700-8700 




22, 24 


2008 June 23-25 


2.2m/CAFOS 


R100 


6000-9000 


2 


21, 23, 24, 25, 29 


2008 December 30 


Palomar/DBSP 


R316 


6200-8700 


10 


4, 6,8 


2009 July 21 


Palomar/DBSP 


R316 


6200-8700 




9, 10 


2009 October 30 


WHT+ISIS 


R316 


5500-8300 




10, 13 


2009 August 03 


CTIO 1.5m/R-C Spec 


47/Ib 


5630-6930 




18,19 



1 ()=low S/N 
2 fiber size 
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Table 4. Stellar parameters for the stars in the sample. 



No. 


Object Id. 


SpT 


Ref 


Ay 


T„ff 


L* 


Mass 


EW[Hck] * 


tt r 1 r\Cr? T 

Ha[10%J 


log M aac 










(mag) 


( K ) 


( L o) 






(km/s) 


( M o/y r ) 


2 


ASR. 118 


K4 


1 


13.0 


4590 


1.17 


1.36 


26 


203 


-10.2 


5 


SSTc2d J032929. 3+311835 


M0 


1 




3850 


0.09 




4.8 






10 


SSTc2d J034227. 1+314433 


K7 


1 


5.0 


4060 


0.22 


0.82 


4.3 






11 


CI* IC 348 LRL 190 


M3.75 


6 


6.5 


3306 


0.09 


0.35 


5 






12 


SSTc2d J034434. 8+315655 


M3 


1 


3.0 


3415 


0.12 


0.45 


130 


504 


-8.0 


17 


Sz 84 


M5.5 


3 


0.5 


3057 


0.15 


0.14 


44 






18 


SSTc2d J161029. 6-392215 


M5 


1 


1.0 


3125 


0.08 


0.21 


18 






19 


RX J1615. 3-3255 


K4 


1 


1.0 


4590 


0.85 


1.28 


20 






22 


SSTc2d ,1162245.4-243124 


M2 


1 


2.0 


3580 


0.29 


0.59 


5.0 


224 




24 


SSTc2d J162506. 9-235050 


M0. 5 


1 


3.0 


3778 


0.21 


0.77 


4.6 


505 


-8.0 


26 


SSTc2d J182829. 1+002756 


K7 


2 


13.0 


4060 


0.34 


0.94 








29 


SSTc2d J182911. 5+002039 


M0 


13 


6.0 


3850 


0.13 


0.72 


25 


050 


-6.5 


.32 


SSTc2d J182935. 6+003504 


K7 


13 


4.0 


4060 


0.52 


1.06 


10.9 


273 


-10.2 


33 


SSTc2d J182936. 2+004217 


F!) 


13 


11.0 


6115 


1.89 


1.30 








34 


SSTc2d J182944. 1+003356 


M0 


13 


5.0 


3850 


0.54 


0.93 


37 


444 


-8.6 


1 


SSTc2d J032835.1+302010 1 " 


K4 


1 


0.0 


4590 


0.07* 




8.1 


449 


-8.5 


3 


SSTc2d .1032903.9+305630+ 


G7 


1 


6.0 


5630 


0.05* 




54 


344 


-9.6 


4 


SSTc2d J032924.1+311958t 


K5 


1 


5.0 


4350 


0.41 


0.96 


6.5 


202 


-10.4 


6 


SSTc2d .1033027.2+302830*'" 


c b 


1 










106 






7 


SSTc2d .1033035.5+311559+ 


K5 


1 




4350 


0.07* 




15.1 


405 


-0.0 


8 


SSTc2d J033038.2+303212t 


M0 C 


1 




3850 


0.05* 






194 




9 


LZK 21 


M0 


1 


6.0 


3850 


1.46 


1.00 


72 






13 


CI* IC 348 LRL 265 


M3 


1 


17. 


3415 


0.42 


0.52 


113 






14 


CI* IC 348 LRL 110 


M0. 5 


1 


5.0 


3778 


0.22 


0.78 


36 


424 


-8.8 


15 


LkHa 329 


K5 


1 


5.0 


4350 


2.89 


1.63 


19.4 


700 


-0.1 


16 


Hn 24 


M0. 5 


1 


2.0 


3778 


0.75 


0.92 


0.2 


311 




21) 


SSTc2d J162148. 5-234027 


M5.5 


5 


7.0 


3057 


0.04 


0.14 


d 






21 


SSTc2d J162221.1-230403t 


K5 


1 




4350 


0.02* 




45.8 


51)4 


-7.1 


23 


SSTc2d J162332. 9-225847 


M0. 5 


1 


5.0 


3778 


0.15 


0.73 


2.8 


413 




28 


DoAr 21 


Kl 


4 


6.0 


5080 


9.84 


1.80 


1-2 


450 




27 


SSTc2d J182858. 1+001724 


G3 


13 


6.0 


5830 


4.76 


1.4 


-2.2 






28 


SSTc2d J182907. 0+003838 


K7 


2 


16.0 


4060 


0.6 


1.04 








30 


SSTc2d J182915. 6+003912 


K5 


13 


12.0 


4350 


0.45 


1.00 


22 


532 


-7.7 


31 


SSTc2d .1182915.6+003923+ 


K7 


2 


0.0 


4060 


0.20 










35 


SSTc2d J182947. 3+003223 


M0 


13 


7.0 


3850 


0.28 


0.85 


49 


618 


-6.9 


Note. 


— For clarity, cold disks arc shown above the horizontal 


line while other objects arc 


below. 







Note. — We used the c2d distances to the clouds to compute the stellar luminosities: 250 pc for Perseus, 187 pc for Cham. II, 150 
pc for Lupus II and III, 120 pc for Ophiuchus and 260 pc for Serpens. 

*Hq equivalent widths arc listed as positive when seen in emission and negative when seen in absorption 

' : Objects identified as likely edge-on disks. Luminosities are presented to show that the stars are underluminous. In these cases, we 
do not calculate masses and ages. 

a : Source is spatially extended in slit. 

k; Continuum spectral type. The object is heavily veiled and no SpT was obtainable at low spectral resolution. 
c : Large uncertainty due to poor S/N and large veiling. 

d : Accreting based on presence of Paschen j3 emission (Gully-Santiago ct al. in prep.). 



References. — 1) Spectral type from this work with optical spectroscopy; 2) Spectral type from this work with photometry; 3) 
iKrautter et al.| i 19971 ); 4) I Jensen et al.| j 20091) 5) Gully-Santiago ct al. in prep. 6) iLuhman et al.| j 20031 ) 13) lQliveira et al.| <2009D 
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Table 5. Optical and near-IR photometry for the cold disk candidate sample 



IN O. 


U E> 


1 7 
V 




ic 


J2MASS 


H2MASS 


K2MASS 




(mag) (mag) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


i 










15.66±0.05 


14.41±0.05 


13.69±0.05 


2 










13.76±0.05 


11.92±0.05 


11.08±0.05 


3 










15.64±0.05 


14.64±0.05 


14.37±0.05 


4 










15.52±0.05 


14.25±0.05 


13.57±0.05 


5 










12.59±0.05 


11.37±0.05 


10.96±0.05 


6 










14.86±0.05 


13.16±0.05 


11.90±0.05 


7 

8 










15.66±0.09 
15.57±0.05 


14.40±0.06 
14.65±0.05 


13.69±0.05 
14.04±0.05 


9 










11.29±0.05 


10.05±0.05 


9.41±0.05 


10 








16.33±0.01 


13.38±0.05 


11.86±0.05 


11.15±0.05 


11 






20.45±0.05 


17.93±0.05 


14.38±0.05 


12.87±0.05 


12.11±0.05 


12 










13.02±0.05 


12.09±0.05 


11.63±0.05 


13 










15.64±0.05 


13.37±0.05 


12.00±0.05 


14 








15.78±0.05 


13.07±0.05 


11.85±0.05 


11.28±0.05 


15 










10.46±0.05 


9.46±0.05 


8.93±0.05 


16 










10.15±0.05 


9.33±0.05 


8.88±0.05 


17 


17.54±0.03 


16.06±0.11 


14.53±0.04 


12.94±0.02 


10.93±0.02 


10.20±0.02 


9.85±0.02 


18 


17.51±0.03 


16.31±0.01 


15.79±0.05 


13.90±0.05 


11.94±0.05 


11.27±0.05 


10.91±0.05 


19 




12.04±0.01 


11.28±0.05 


10.54±0.05 


9.55±0.05 


8.78±0.05 


8.56±0.05 


20 






19.59±0.04 




13.59±0.05 


12.35±0.05 


11.69±0.05 


21 


19.18±0.03 


17.74±0.01 


17.03±0.04 




14.66±0.05 


13.87±0.05 


13.57±0.05 


22 




15.74±0.50 


14.82±0.05 


12.51±0.02 


10.37±0.05 


9.43±0.05 


9.06±0.05 


23 






16.22±0.04 




11.49±0.05 


10.68±0.05 


10.29±0.05 


24 


18.38±0.05 


16.59±0.05 


15.55±0.05 




11.06±0.05 


10.03±0.05 


9.51±0.05 


25 


15.80±0.75 


14.00±0.50 


12.00±0.05 


10.05±0.90 


8.09±0.02 


6.86±0.05 


6.23±0.02 


26 










14.99±0.05 


12.76±0.05 


11.81±0.05 


27 










10.72±0.05 


9.83±0.05 


9.45±0.05 


28 










15.20±0.05 


12.90±0.05 


11.55±0.05 


29 










14.03±0.05 


12.89±0.05 


12.42±0.05 


30 










14.53±0.05 


12.65±0.05 


11.69±0.05 


31 






18.38±0.05 




14.14±0.05 


12.79±0.05 


12.13±0.05 


32 










12.00±0.05 


10.84±0.05 


10.26±0.05 


33 










13.27±0.05 


11.92±0.05 


11.31±0.05 


34 






19.04±0.50 




12.19±0.05 


11.15±0.05 


10.77±0.05 


35 






19.14±0.05 




13.45±0.05 


12.11±0.05 


11.45±0.05 



References. — Optical fluxes come from 1) Bouy et al. (in prep) f or Pe rseus sources (1-15); 2) ISpezzi etHl j2008h 
for the only object from Chamaeleon II (#16); 3) IComeron et a hl d200 9T) f or the Lupus sources (# 17 and 18); 4) 
Bouy et al. (in prep.) for Ophiuchus sources (19 - 251: and 5l lSpezzi et all l|2010l) for the sources in Serpens (26-35) 
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Table 6. Selection criteria for edge-on disks 



No. 


Under 


Forbidden 


Silicate 


Ice 


High 




luminous 


lines 


absorption 


absorption 


L Disk/ L * 


1 










✓ 


3 


✓ 


■n 


s 


s 




4 






■/ 


s 




6 












7 


✓ 




s 




■f 


8 


✓ 




s 




■r 


21 


✓ 


s 






■r 


31 










/? 
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Table 7. Disk physical parameters and observables 



No. 


Reg. 


(AU) 


Inc. 

(deg) 


Mass 1 " 
10" 3 Mq 


Sett. 


F mm 

(mJy) 


M disk 
10-3 M 


Ref. 


F30/F13 


■Ldisk/L* 


nRob. 
re. 1 

hole 

(AU) 


Model ID 
(Rob.) 


1 




edge 








<1.0 


<0.3 


1 


20.17 








2 


B 


1.5±1 


75 


obs 


0.5 


3.7±0.9 


1.1±0.3 


1 


4.29 


0.16 


1.64 


3005228 


3 


B 


edge 














21.35 








4 




edge 








<1.3 


<0.4 


1 


9.82 








5 


A 


g+4 
°-2 


10 


obs 


0.5 


6.3±1.1 


1.9±0.3 


1 


12.86 


0.2 


15.3 


3000829 


6 




edge 








<0.7 


<0.2 


1 


17.63 








7 


B 


edge 








1.3±0.5 


0.4±0.2 


1 


17.47 








8 




edge 














33.19 








9 


B 


no hole 


30 


obs 


0.5 


13. Oil. 3 


4.0±0.4 


1 


3.01 


0.19 






10 


A 


5±1 


60 


0.001 


1 


<1.4 


<0.4 


1 


4.76 


0.04 


1.87 


3018041 


11 


A 


5±1 


30 


0.003 


1 








7.85 


0.22 


0.75 


3015642 


13 




no hole 


15 


obs 


1 


11.6±1.0 


3.5±0.3 


1 


4.51 


0.18 


1.71 


3007226 


12 


A 


3±1 


5 


0.002 


1 


<1.4 


<0.4 


1 


5.0 


0.2 


1.04 


3017747 


14 


B 


no hole 


60 


0.3 


1 


<1.5 


<0.5 


1 


5.78 


0.23 


0.8 


3019220 


15 


B 


no hole 


75 


0.005 


1 


<1.2 


<0.4 


1 


2.38 


0.05 


8.6 


3005263 


16 


B 


no hole 


65 


0.05 


1 








2.65 


0.12 




3012969 


17 




55±5 


45 


0.03 


1 


<35.9 


<4.0 


2 


11.47 


0.15 


96.0 


3004704 


18 


A 


2±1 


75 


0.5 


0.05 








4.67 


0.09 




3000691 


19 


A 


2±1 


5 


0.5 


0.1 








5.36 


0.09 


3.8 


3017818 


20 




no hole 


5 


0.01 


1 


<4.6 


<0.3 


1 


3.82 


0.33 


0.8 


3019048 


21 




edge 








<4.4 


<0.3 


1 


5.55 








22 


A 


L -0.5 


30 


0.25 


1 


<5.4 


<0.4 


1 


3.16 


0.18 


2.75 


3010142 


23 


B 


no hole 


80 


0.004 


1 








3.55 


0.19 


0.3 


3004623 


21 


A 


3±2 


5 


0.1 


0.25 








6.98 


0.13 


2.21 


3005973 


25 


A 


back 














5.67 


0.12 






26 


A 


4±1 


5 


obs 


0.15 


24.0±0.6 


8.0±0.2 




8.87 


0.17 


17.6 


3002156 


27 




back 














1.28 


0.03 






28 




no hole 


5 


0.002 


1 








2.41 


0.15 




3008873 


20 


A 


8±2 


30 


0.01 


1 








11.71 


0.07 


3.76 


3010258 


30 




no hole 


15 


0.05 


1 








2.53 


0.03 


0.5 


3006320 


31 




edge 


81.4 












5.89 








32 


A 


7±1 


5 


0.01 


1 








6.01 


0.11 


4.14 


3008745 


33 


A 


13±3 


5 


obs 


0.005 


2.3±0.6 


0.8±0.2 


1 


7.36 


0.33 


3.94 


3018536 


34 


A 


25^ 


5 


obs 


0.25 


6.3±0.6 


2.1±0.2 


1 


10.63 


0.02 


8.94 


3019329 


35 




no hole 


80 


obs 


0.025 


2.2±0.9 


0.7±0.3 


1 


1.5 


0.13 


0.4 


3016446 



References. — 1) C onti nuum flux at 1 . 3 mm fro m this work fij 12.41: 2) iNuernberger et al ] (jl997h : 3) lOsterloh fc Beckwithl l |l995h : 
4) lAndre et al] <1990h : 5) iHenning et all Jl993h : 6) ISvlvester et all jl996h . 

Note. — Column "Ref." gives the regions in the upper left panel of Figure [TS] 

Note. — * : "edge" labels edge-on disks; "back" labels the two SEDs contaminated by background cloud material. 
Note. — t; obs indicates that the disk mass in the model is the observed millimeter mass. 
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Table 8: Classification and mineralogy of the disk sample 



No. 


Name 


Hole 


Acc. 


Mdisk 


10 flM 


20 fim 


PAH 


Cryst. 


2 


ASR118 


Y 


Y 


L 


Y 


Y 


N 


T 


5 


SSTc2d J032929. 3+3118 


Y 


N 


L 


Y 


Y 


N 


T 


10 


SSTc2d J034227. 1+3144 


Y 


N 


L 


Y 


N 


N 


Y 


11 


CI* IC 348 LRL 190 


Y 


N 




Y 


Y 


Y 


N 


12 


SSTc2d J034434. 8+3156 


Y 


Y 


L 


Y 


Y 


N 


Y 


17 


Sz 84 


Y 1 


Y 


L 


N 


N 


N 


N 


18 


SSTc2d J161029. 6-3922 


Y 






Y 


Y 


N 


Y 


19 


RX J1615.3-3255 


Y 


Y 




Y 


Y 


Y 


N 


22 


SSTc2d J162245.4-2431 


Y 


N 


L 


Y 


Y 


N 


N 


21 


SSTc2d J162506. 9-2350 


Y 


Y 




Y 


Y 


N 


Y 


26 


SSTc2d J182829.1+0027 


Y 




11 


Y 


Y 


N 


T 


29 


SSTc2d J182911. 5+0020 


Y 


Y 




Y 


Y 


N 


N 


32 


SSTc2d J182935. 6+0035 


Y 


Y 




Y 


Y 


N 


Y 


33 


SSTc2d J182936.2+0042 


Y 




L 


Y 


Y 


N 


T 


34 


SSTc2d J182944. 1+0033 


Y 


Y 


L 


Y 


N 


N 


N 


1 


SSTc2d J032835.0+3020 1 " 


N 


Y 


L 


Y 


Y 


N 


N 


3 


SSTc2d J032903.9+3056t 


N 


Y 




N 


Y 


Y 


Y 


4 


SSTc2d J032924.1+3119t 


N 


Y 


L 


N 


Y 


N 


N 


6 


SSTc2d J033027.1+3028t 


N 


Y 


L 


Y 


Y 


N 


Y 


7 


SSTc2d J033035.5+3115t 


N 


Y 


L 


N 


Y 


N 


Y 


8 


SSTc2d J033038.2+3032t 


N 


Y 




K 


Y 


N 


Y 


9 


LZK 21 


N 


Y 


L 


Y 


Y 


N 


Y 


13 


CI* IC 348 LRL 265 


N 




L 


Y 


Y 


N 


T 


11 


CI* IC 348 LRL 110 


N 


Y 


L 


Y 


N 


Y 


N 


15 


LkHa 329 


N 


Y 


L 


Y 


Y 


N 


Y 


16 


Hn 24 


N 


N 




Y 


Y 


N 


Y 


20 


SSTc2d J162148.5-2340 


N 


Y 


L 


Y 


Y 


N 


N 


21 


SSTc2d J162221.1-2304t 


N 


Y 


L 


Y 


Y 


N 


N 


2.3 


SSTc2d J162332.8-2258 


N 


N 




Y 


N 


N 


Y 


25 


DoAr 21 


N 


N 


L 


N 


N 


Y 


N 


27 


SSTc2d J182858.1+0017 


N 


N 




N 


N 


Y 


N 


28 


SSTc2d J182907.0+0038 


N 






Y 


N 


Y 


N 


30 


SSTc2d J182915.6+0039 


N 


Y 




Y 


Y 


N 


Y 


31 


SSTc2d J182915.6+0039t 


N 






Y 


Y 


N 


N 


35 


SSTc2d J182947.3+0032 


N 


Y 


L 


Y 


Y 


N 


N 



Note. — For clarity, cold disks and other objects are shown above and below the horizontal line, respec- 
tively. Y indicates feature is present, N indicates that the feature is not present, and T indicates a tentative 
detection. 

Note. — 1 : The classification of Sz 84 (# 17) as a cold disk or an extended source is ambiguous but we 
keep it in assuming it might be both. 



Note. — t; Close to edge-on according to SED fit and luminosity below Main Sequence. In these cases, 
we do not classify as cold disks, since the hole size determination from the SED fit is highly uncertain. 
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5 10 15 20 25 30 35 
A (jam) 

Fig. 1. — IRS spectra of the of the close to edge-on systems ordered from bottom to top with 
growing effective temperature. The thick and grey lines are the binned and original spectra, 
respectively and the numbers give their identifications in Table [TJ 
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Fig. 2. — IRS spectra of the disks in the sample for which inner holes were not required to 
obtain a good SED fit organized from bottom to top with increasing SED slope a. The thick 
and grey lines are the binned and original spectra, respectively and the numbers give their 
identifications in Table [TJ 
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5 10 15 20 25 30 35 
A (fxm.) 

Fig. 3. — IRS spectra of the disks with potentially extended mid-IR emission from surround- 
ing clouds ordered from bottom to top with growing stellar mass. Object 17 (Sz 84) is shown 
here as comparison due to its similar IRS spectrum, although it is classified as a cold disk. 
The thick and grey lines are the binned and original spectra, respectively and the numbers 
give their identifications in Table [TJ 



5 10 15 20 25 30 35 
A (yu.m) 

Fig. 4. — IRS spectra of the accreting cold disks in the sample ordered from bottom to top 
with growing inner hole radii. The thick and grey lines are the binned and original spectra, 
respectively and the numbers give their identifications in Table [TJ 
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i i i i I i i i i I i i i i I i i i i I i i i i I i i i i I i i i i 

5 10 15 20 25 30 35 
A (yu.m) 

Fig. 5. — IRS spectra of the non-accreting or non-classified cold disks ordered from bottom 
to top with growing inner hole radii. The thick and grey lines are the binned and original 
spectra, respectively and the numbers give their identifications in Table HJ 
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Fig. 6. — Spectral Energy Distributions of the eight close to edge-on disks (top) in the sample 
and the two extended sources (bottom). Points are observed fluxes from IRAC, MIPS and 
IRAM. The thick solid line is the IRS spectrum and the numbers give their identifications 
in Table [TJ 
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Fig. 7. — Spectral Energy Distributions of the disks with holes in the sample. Points are 
dereddened fluxes from IRAC, MIPS and IRAM. The thick solid line is the IRS spectrum, 
the dotted line is the stellar photospheric model, the red dash-dot line is the best-fit RADMC 
disk model, and the thin black dashed line is the Robitaille model fit. The numbers give 
their identifications in Table [TJ 



-51 - 



E 



-9 - 

-10 - 

-11 -$ 
f 

-12 I 
-13 - 
-14 - 
-15 - 
-9 - 
-10 - 

. 

-12 f 
-13 - 
-14 j- 
-15 - 
-9 - 
-10 - 

-12 f 
-13 - 
-14 - 
-15 - 
-9 - 
-10 - 
-1 1 
-12 
-13 
-14 
-15 



22 



4 



\ T 



26 



-4 

_L 



32 



34 



24 



29 



33 



1 10 100 

Wavelength (/xm) 



1000 



1 10 100 1000 

Wavelength (/Am) 



Fig. 8. — Figure 7 continued. 
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Fig. 9. — Spectral Energy Distributions of the disks without holes in the sample. Points are 
dereddened fluxes from IRAC, MIPS and IRAM. The thick solid line is the IRS spectrum, 
the dotted line is the stellar photospheric model and the red dash-dot line is the best-fit 
RADMC disk model. The numbers give their identifications in Table [TJ 
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Comparison b etween hole radii computed with the RADMC models 
( Dullemond fc Dominik 2004 ) and with the iRobitaille et~aL ( 2007 ) disk model grid. Even 



Fig 



10.- 



with the large scatter, both quantities are correlated, which shows that the derived disk hole 
radius is model independent to generally within a factor of 2-3. 
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Fig. 11. — Inner disk radius histograms for the cold disks in the sample. The cold disks 
predominantly have inner holes with radii smaller than 10 AU. Many of the cold disks in the 
literature have larger hole sizes. 
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Fig. 12. — Continuum subtracted PAH emission from object # 27. The vertical dashed lines 
show PAH features. 



-56 - 




Fig. 13. — Shape-strength diagram for the 10 /iin silicate emission feature for all cold disks in 
this s ample (filled dots), compared with those of the c2d IRS T Tauri sample (jOlofsson et al. 
20091 . grey dots). 
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Fig. 14. — Histog ram of F30/F n rat ios for the disks in this paper, with the definition of 
the boundaries by iBrown et al.l ( 120071 ) marked by vertical dashed lines for comparison. The 
region on the left is supposed to contain disks without holes, the middle region holds cold 
disks and on the right edge-on disks. The actual distribution of these types of sources in 
our sample is shown with normal disks in the blue short dashed line, cold disks with solid 
black and edge-on disks with the red long dashed line. There is significant overlap between 
all categories around F 30 /F 13 of 5 but the general trends hold true in our sample. 
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Fig. 15. — Selection criteria. Red stars and crosses are the spectroscopically confirmed 
cold disks, black diamonds are disks without holes and bl ue triangles are ed ge-on disks. 
The grey crosses are all the sources in c2d YSO catalog (jEvans et all 120091 ) . (top left) 
Our selection criteria. Regions A and B define the two selection criteria for identifying 
cold disks in photometric sample, the first one, represented with red crosses, selects "clean" 
inner holes and the second one, represent ed with red s tars, identifies cold disks with near-IR 
exces s, (top right) Se l ection criteria from Fang et al.l ( 120091 ) . (botto m left) Selec t ion cr iteria 
from iMuzerolle et al.l (120101 ) . (bottom right) Selection criteria from ICieza et al.l ( 120101 ) with 
transitional disks according to their definitions in the upper left quadrant. Arrows show the 
effect of correcting for an extinction of Av=10. 
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Fig. 16. — Correlation between hole size and disk mass. Dots and upper limits represent 
source with millimeter fluxes from t his work while th e dia monds are the cold disks with mil- 
limeter measured disk masses from iKim et al.l (120091) andlBrown et al.l (120071 ) . To minimize 
systematic differences, disk masses for sources fromlKim et al.l (120091 ) have been recalculated 
using the 1.3 mm fluxes in lAndrews &: Williams! ( 120051 ) and the conversion in Section 14.11 
The correlation has a Pearson's correlation coefficient of 0.6 and is statistically significant 
on the 99% level including the upper limits and 90% level with only the detected sources. 
The systematically lower disk masses in this survey might explain the generally smaller hole 
sizes than previously published transitional disks. 
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Fig. 17. — Relationship between the size of the inner disk hole and the mass of the central 
star. Solid circles repre sent our sample while stars and diamonds are transitional disks 
from iBrown et al.l (120071 ) and Kim et al.l ( 120091 ) respectively. The dotted line is the best 
linear fit to the data. The small crosses mark the semimajor axis of all exoplanets in the 



http://exopl.anet .eu database as of December 2009, compared to the mass of their parent 



stars. 
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Fig. 18. — Relationship between the size of the inner disk hole and the intensity of the 10 
/jm feature. Solid, grey and open circles represent c2d accreting, non-accreting cold disks 
and objects for which we do not have yet accretion information, respectively. 



